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In this work, we study the Variational Iteration Method (VIM) accompanied by He's polynomials for obtaining the 


solution of the temperature distribution in the casting-mould heterogeneous system with fractional order. The combined 
form of VIM and the Homotopy Perturbation Method (HPM) is called the Variational Iteration Homotopy Perturbation 
Method (VIHPM) which yields to implement the system of equations directly. The identification of the Lagrange 
multiplier is essentially more reliable and higher accurate for such type of problems. The approximate solution converges 


rapidly to the exact solution which confirms the accuracy of this approach. Some graphical representations are 


demonstrated to show the validity of this approach. 
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Introduction 


Partial differential equations play an important and expanded significance in describing a variety of 
physical phenomena in natural phenomena, such as solid-state physics, fluid mechanics, plasma 
physics, population dynamics, chemical kinetics, nonlinear optics, protein chemistry, and so on. For 
appropriate topics, these nonlinear models, as well as their analytic solutions, are quite interesting. 


He [1], [2] developed semi-analytical methods such as the Variational Iteration Method (VIM) and 
the Homotopy Perturbation Method (HPM) for the first time to obtain the approximate solution to 
linear and nonlinear problems. Later, many scientists confirmed that HPM is a very powerful 
technique as Edalatpanah and Rashidi [3] confirmed the validity of HPM for finding the approximate 
solution of the convection-diffusion equation. Yildirim and Ozis [4] studied HPM for the 
approximate solution of the Lane-Emden type problem. Domairry and Aziz [5] used HPM to obtain 
the approximate solution of MHD squeeze flow between two parallel disks. We recommend the 
researchers study the efficiency of HPM in [6]. 
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Grzymkowski et al. [7], [8] employed HPM whereas Tripathi and Mishra [9] adopted HPM together with 
the Laplace transform to determine the temperature distribution in the casting-mould heterogeneous 
system as a continuous function, which is particularly useful for analyzing the mould. Vanani et al. [10] 
used a weighted approach based on HPM to solve the heat equation in the cast-mould heterogeneous 
domain. Later, this proposed approach has also been examined in more than one spatial dimension, 
indicating that this method has a broader application in nonlinear PDE systems [11], [12]. This study is 
particularly powerful for fractal theory and fractal calculus, and it can be seen as dependable in getting 
analytical solutions and suitable for other nonlinear issues [13]-[15]. 


This study presents the idea of Variational Iteration Homotopy Perturbation Method (VIHPM) to obtain 
the solution of casting-mould heterogeneous system. The quality of the current method is appropriate to 
provide the analytical results to the given examples. This study is summarized as follows: In Section 2, we 
construct the description of the casting-mould system. In Section 3 and 4, we examine the idea of VIM 
and HPM. We formulate the idea of VIHPM for the implementation of casting-mould system in Section 
5. In Section 6, we perform this scheme to a numerical problem to show its capability and efficiency 
whereas the results and discussion along conclusion are discussed in Section 7 and Section 8. 


2 | The Description of the Problem 


In this frame, we develop a casting-mould system to examine temperature distribution. Consider, O(¢, č) 
for casting and 3(c,€) for mould be two regions on the boundary of the problem [16]. 


(S) = (c, E): G € EP 0], 2 € [0, e (1) 
= =(c,&):¢ € [0,cal, & € [0, €). (2) 
These two functions O(Ç, €) and 3(C, &) satisfy the heat transfer equation inside the domains such as 


IOE) _ 9° O(c, 8) 
a6 A Jg 


, (é) E0, 


(3) 
aze) _ 07K (c, &) 
d 7? Jg 


Where a, and a, are the thermal diffusivity, © and £ represent the temperature, ç and € show the spatial 


; (GE) EB: 


time respectively with the following conditions 
O(c,0) = =1(c), on Ty 
B(c,0) = E2(c), onT'4 
O(c, 6) =O(E),onT3 


Se = a(),onTs ) 
©(0,6) = (0,6), onT2 

00, 20,6 
ò AA = ð onl’? 


The determination of the approximate solution of the casting-mould system requires careful selection of 
these boundary parameters. 


3 | Variational Iteration Method (VIM) 


The VIM was first proposed by He [1], [16] and Nadeem and Li [17]. This method is now widely used in 
many fields such as Physics [18], Chemistry [19], Biomedical [20] and engineering sciences to study linear 
and non-linear partial differential equations [21]. Consider a differential problem such as 


LO(c, €) + NO(S, €) = g(<, §). (5) 


ok. ; 
Where L = zp ÍS linear and N be a nonlinear operator and g(c,é) represents the source term. The VIM 


examines the correction functional as follows, 


Onl, E) = Onlc, E) + I “ Ms)[LO,(C,s) +NOq(<,s) - g(s) as. 6) 


Eq. (6) is called a correction functional, @p is called initial approximation, A denotes Lagrange multiplier 


is obtained using variational theory and O, is called restricted variable such that 60, = 0, so 


A= -1, fork = 1, 
A=s—E fork =2, 


And in general, for m > 1 


_ He- i 


DE T) 
Eqs. (6) and (7) develop the following iteration formula, 
E (-1)™(5 — k-1 = 
BrE Da f° TS e + NSl) -ges 


Thus, Eq. (8) is said to be a correction functional formula. The values of the Lagrange multiplier A will 
be identified optimally via integration by parts. The consecutive approximation ©,,41,n 2 0 of the 
solution © can acquire with the help of the Lagrange multiplier. Ultimately, the solution is given by O = 
limpo Os. 


4 | The Basic Idea of the HPM 


Let us consider the following non-linear differential problem [22], [23]: 


T(©) — g(w) =0. 9) 
With boundary conditions 
(0, £) =0. (10) 


Where T is a functional and t is a boundary operator, 9(w) is a known function. We separate T into two 
operators, such as T; is a linear and T, represents a nonlinear operator respectively. Thus, Eg. (9) 
becomes as 


L(©) + N(©) - g(w) = 0. (11) 
Let's consider a homotopy £(r, p):Q x [0,1] > R that satisfies 

H(E,p) = (1 - PILE) - L(@p)] + pIL() - N(E) - g(w)]. (12) 
Or 

H(2,p) = L(€) - L(©o) + pL(©o) + pIN() - g(w)] = 0. (13) 


Here p € [0,1], is an embedding parameter, and Op is an initial guess of Eg. (9), which satisfies the 
boundary conditions. On the other hand, HPM considered p as a small parameter, and suppose that the 
solution of Eg. (13) can be termed as a power series 


E = Ho + pv +p*Ey +o. (14) 


Put Eg. (14) into Eq. (13) and comparing the identical powers of p, we can find a series of equations in 


such a form 


po E -H(So, 24 =0 
p®:23-H(So,21,82) = 


Where H(So, By, Ezg, z) depend upon Eo, 21, E2, +: poy and subjected to He's polynomials. We may 


calculate them by using the formula 
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H (Bo Er 820,5) = 5N (Zico Sip')| (16) 


Let p = 1, the analytical solution of Eg. (9) is as follows 


@=limS = 89+ 2, +8, + 834°. (17) 
pol 


It is obvious that the system of nonlinear equations in Eg. (15) is easy to solve, and the components &;,i > 
0 of the HPM can be completely determined, and the series solutions are thus entirely determined. 


5 | VIM with He's Polynomials (VIHPM) 


In this segment, we construct the strategy of the VIMHP. Therefore, consider the following equation 


L, O(c, é) + N, O(c, é) = g1(c, €) 


= = ; 18 
L2B(c, €) + N2B(c, €) = 82(6, 8) a 
According to VIM, we can write the correction functional as follow 
€ 
Ons (6E) = Oale, č) + [7 Ar()[Liun(c,s) + Ni@n(c,s) - g1(¢,s)]ds m 


[m 


n1 (CE) = Eal E) + [° ho(s)[LiEa(G,8) + NoEa(s,8) - g2(678)]ds 


Here A, and A, are Lagrange multipliers. Now, using HPM [24], [25], a system of Eg. (19) becomes as 
(oe) n = & co n 
Zo P Onl, E) = O0, E) + p [7 MONEZ p'O) - g1(c,8)lds 
(oe) n EEE € co no : 
n-0 P"On(G €) = Eo(<, €) + pi Ag(s)IN2(n=0 P"=En(<,S)) — S2(c,s)]ds 


Equating the components of Eg. (20) with similar powers of p, and using p approaches to 1, we can obtain 


(20) 


following the formula. 


O(c, é) = Oo(c, é) + Ox(c, č) + O2(c, é) pees limp 1 Yn-1P"On(s, &). (21) 
Be é) = Bole; é) + Ox(c, é) + Balc, é) tre limp +1 bye en (ee é). (22) 


6 | Numerical Applications 


In this segment, we implement the idea of VIHPM for the approximate solution of casting-mould 
heterogeneous system. The significant findings demonstrate the validity and consistency of this scheme. 


6.1 | Test Problem I 


Considering a, = 7 a =1, 6, =1, ô, = 2, we get the system of Eg. (3) such as, 
d© _10°@ 08 078 (23) 
IE 40c2 7 IE Ac?” 


With initial condition 
OQ(c,0) =e*5, E(c,0) =e°. (24) 


According to the strategy, the correctional functional for the system of Eg. (23) becomes as 


Orale E) = One 8) + S mo e- Seas. 25) 


B 


Ens (cr 6) = Enlc 8) + f M9) P- Sas 


Which implies the stationary terms 


1+A,(s)=0, Aj(s=é) =0, 
1+A ,(s)=0, Ag(s=€) =0 


The Lagrange multipliers yields 
Aa(s) = Ag(s) = - 


Put these values into the system of Eg. (25), 


1320, 
On (GE) = O(G E) - Je Se l 


n dEn 
Eniil(G E) = EnG €) -f É ðs = = 


Applying the VIHPM on the system of Eg. (26), we get 


(26) 


co 


Y Pre = O68) +P f k 2 p” a ds, 


n=0 


> Ps pre BCDE +P f ae 


Now from Eq. (23), we can select @o(C, €) = e°, Zy(C, &) = eS. Accordingly, the coefficients of p are, For 


zero-th order of p, we get 


p°: Oolc, č) = e% 
p°: EAC é) = et l 


For the 1st order of p, we get 


1940 
pata: f f as Sate 


0-8 
pt: £) = f Eees = ere 
0 


For the 2.4 order of p, we get 


1020 
p?:0,(¢,&) = Cs 388 =5 -e22 


94E 
p°: E2(c,&) = ce Jase! zee 


For the 3*1 order of p, we get 


t [1320 1 
poa f f las = eri 


4 dG 6 
t [028 1 
3-7 = =2 — Leg 


For the 4 order of p, we get 
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E 2 
ptos f R 3 as = ere! 


4 dc2 24 
č eka) 1 l 
4.7 = 3 = cr4 
pP 24c, $) J E Jos oT ag a 


The series solutions are therefore given by 


E2 E3 E4 
OGE = Op +; +O, 4054 =e (IHEH TH Tea], 
E M z _ E E2 E3 z1 
B(s, E) = Eo +81 +23 +E3 ++ we'(ises eiit í 


And hence the exact solution becomes 


O(Gs) =e", Mose, 


6.2 | Test Problem II 


Next, considering the system of Eg. (22) with the following initial conditions 


O(c, 0) =2+e%, E(c,0) = eS. (27) 


The identification of Lagrange multipliers will be the same and thus according to VIHPM strategy, the 
coefficients of p are, For zero-th order of p, we get 
p?:Oo(C,) =2+e%, 
p°: Eol, 6) = e. 


For the 1* order of p, we get 


(28) 
£ [1920 
pots f f aaas = 2+ exe 


£ [o2z 
paas f elas = ect 
0 


For the 2" order of p, we get 


č [1920 1 
p?:@(C, €) = J f z s= 


£928 1 
pE: &) = f | zilas = ~e&? 
0 


For the 34 order of p, we get 


£17122 
p*: Oa(c, €) = f f : Sas = Lag 
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4 0c2 6 
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For the 4th order of p, we get 


40 

é [925 1 
4.0 = 3 eG 
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The series solutions are therefore given by 


ae ee 
O(c, €) = On + ©, +0: + O34 --: -2re* (eer er cece | 
2 3 4 
Els, E) = Ho +21 +82 +23 +- sefiri iriti | 


And hence the exact solution becomes 
T [cE =2+e%", Be, 8) =e, 


| Results and Discussion 


In this segment, we discuss the results of casing-mould system obtained by VIHPM. We see that the 
identification of the Lagrange multiplier in both examples are very easy to find. We investigate the 
solution results in the form of the series which converges rapidly to the exact solution in both examples. 
We observe that the change in the initial condition of casting system does not affect to the mould system 
which show the validity of VIHPM. 


8 | Conclusion 


In this work, we implement a semi-analytical approach VIHPM to find the approximate solution of 
casting-mould heterogeneous system together with intimal conditions. In this work, first, we utilize VIM 
to construct the correct function for the identification of the Lagrange multiplier and then HPM using 
He's polynomials has been introduced to get the iterative formula. This strategy is very helpful for 
finding the series solution of the linear system. It is important to note that the VIHPM has a high rate 
of convergence towards the exact solutions. Eventually, we observed that the suggested approach is very 
strong, competent, and convenient in obtaining analytical solutions. We recommend this strategy for 
several other nonlinear cases, particularly fractal calculus and fractional calculus in science and 
engineering. 
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